Magnetic Phase Diagram of FeAs based superconductors 
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The recently discovered high-temperature superconductivity in doped quaternary iron oxypnic- 
tides correlates experimentally with a magnetic instability. We have used first-principles calculations 
to determine a magnetic phase diagram of ReO i-s FeAs (Re=La-Dy) as a function of the doping 
8, of the FeAs in-plane lattice constant a, and of the distance between the Fe and As planes, that 
is qualitatively consistent with recent experimental findings on the doping, internal (chemical) and 
external pressure dependence. The existence of a tricritical point (TCP) in the phase diagram 
suggests new ways of enhancing T c . 

PACS numbers: 



The recent discovery of superconductivity in electron- 
doped La[Oi_ x FjFeAs with a critical temperature (T c ) 
of 26 K [1 has stimulated a massive experimental 
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to find other, higher T c materials in this completely 
new family of iron-pnictide superconductors. The un- 
doped parent compound LaOFeAs is a poor metal with 
an ordered antiferromagnetic ( AFM ) ground state [5j [11] 
but with increasing F doping [12] the magnetic order- 
ing is suppressed and superconductivity emerges. This 
strongly suggests that magnetic fluctuations in the iron 
layers close to the quantum critical point (QCP) play a 
fundamental role in the superconducting pairing mecha- 
nism. Recent experimental [I2j [I3j [14] [29] and theoret- 
ical [25 , 26 , 2TJ EH] results also suggest that spin fluctu- 
ations in the vicinity of the QCP mediate the supercon- 
ductivity as in the cuprates, heavy fermion materials, or 
ruthenates [30] [31] [32] . In view of this, so far empirical, 
correlation between superconductivity and QCPs, it is 
important to understand how the magnetic (in) stability 
depends on structural and chemical parameters that are 
accessible to experiment. That is the subject of this pa- 
per. 

Superconducting ReOFeAs crystallizes in a tetragonal 
layered structure with PA/nmm symmetry and consists 
of layers of covalently bonded FeAs alternating with lay- 
ers of more ionically bonded ReO. With eight atoms (two 
formula units) in the unit cell, it can be described us- 
ing only two internal structural parameters in addition 
to the lattice parameters a and c 33]. One of the inter- 
nal parameters, dpe-As? describes the separation between 
planes of Fe and As; the other, <iR e _o, between planes 
of Re and O. The excess electron from the ReO layer 
is donated to the FeAs layer which is metallic with a 
number of partly filled bands of mainly Fe 3d character 
intersecting the Fermi energy p~9j [34] . In addition to the 
four structural parameters, we define a doping 5 in terms 
of deviations from this ideal stoichiometry per Fe atom. 
Electron and hole doping can be achieved, for example, 



by partially replacing oxygen with fluorine [T] or trivalent 
La with divalent Sr [35] . 

Since mapping out the magnetic phase diagram in five 
dimensions is impossible, we need to identify a smaller 
number of key independent variables. Clearly the doping 
S is one. Although the interaction between the ReO and 
FeAs layers is by no means negligible, it is widely ac- 
cepted that the superconducting properties of ReOFeAs 
emerge from the FeAs layers. The main role of the ReO 
layers is to determine the lattice parameters [33] and to 
contribute doping electrons while the magnetism and su- 
perconductivity are associated with the "active" FeAs 
layers so we fix c?R e -o and c at their calculated equilib- 
rium values for undoped LaOFeAs. The in-plane lattice 
constant a and the separation between the Fe and As lay- 
ers (ipe-As ED HH US [24] are the other two key parame- 
ters we identify to construct magnetic phase diagrams in 
the (5, a = a eg , <iFe-As) and (S = 0.25, a, dp e -As) planes 
of this parameter space. We do this by calculating the 
energy E(8, a, dpe-As) f° r non-magnetic (NM), checker- 
board AFM (C-AFM) and stripe AFM (S-AFM) order- 
ings, going considerably beyond existing attempts to de- 
termine the QCPs of FeAs-based materials as a function 
of doping or pressure [25] [28] . 

The total energy calculations are carried out within 
the framework of density functional theory (DFT) us- 
ing a spin polarized generalized gradient approximation 
(SGGA) for the exchange-correlation potential (PW91 
functional). The electronic ground state is calculated by 
solving the Kohn-Sham equations self-consistently with 
the projected augmented wave method (PAW) [36] and 
a cut off energy of 500 eV for the plane wave basis as 
implemented in the Vienna ab initio simulation package 
(VASP) [37]. The C-AFM and S-AFM states are de- 
scribed in a y/2 x y/2 x 1 tetragonal structure. Though 
an orthorhombic distortion is observed [5] [21] for the S- 
AFM ordered state in the LaFeAsO parent compounds, it 
is suppressed by doping. Since we are mainly interested 
in locating the phase boundaries, it can be neglected. 
The Brillouin zone integrations are performed with the 
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FIG. 1: Phase diagram of LaOFeAs as a function of dop- 
ing S and the distance between the Fe and As planes along 
the c axis, dpe-As- Blue, red and green represent respectively 
checkerboard AFM (C-AFM), stripe AFM (S-AFM) and non- 
magnetic (NM) ordering. In the AFM region, more stable 
magnetic ordering is described by darker color, do is the cal- 
culated equilibrium value of dFe-As for the undoped material. 
On the right hand side we show the energies per Fe atom of C- 
AFM, S-AFM and NM ordered states relative to their average 
value as a function of doping for the vertical line cfce-As = do. 
The white arrows describe the pressure dependence of over- 
and under-doped LaOFeAs. 



improved tetrahedron method [38 with a sampling grid 
of 6 x 6 x 4 k-points plus T point. When using pseu- 
dopotentials, it has been shown that it is necessary to 
treat the 3p states of Fe as valence states to reproduce 
all-electron results [24]. Doping was modelled by adding 
or substracting electrons and compensating their charges 
with a homogeneous fixed background charge [37]. This 
approximation was checked explicitly using supercell cal- 
culations in which a fraction of the oxygen atoms were 
replaced with fluorine, and proved not to be critical. 

We begin by determining the ground state of LaOFeAs 
as a function of S and dpe-As- First of all the energy is 
minimized with respect to the internal structural parame- 
ters c^Fe-As and d^e-o for the undoped parent compound 
LaOFeAs using the experimental a and c [39] . This opti- 
mized structure is then frozen and the total energy cal- 
culated as a function of dpe-As and 5 for NM, C-AFM 



and S-AFM ordering. The phase diagram obtained from 
these energies is shown in Fig. [I] where positive and neg- 
ative S correspond to electron and hole doping, respec- 
tively. We propose that this phase diagram describes 
qualitatively all FeAs based materials since substituting 
different elements in the Re-0 layers only changes the 
effective doping S and the lattice parameters a and c. 
The effect of changing c, for example by external uniax- 
ial pressure, is given in our phase diagram in terms of its 
effect on dpe-As; changing the Re position has otherwise 
little effect. The effect of changing a will be discussed 
below. Because our phase diagram is calculated as a 
function of two variables, the boundaries separating re- 
gions with different ground states, shown as dashed lines 
in Fig. [I] are in general quantum critical lines (QCL). 
The uncertainty in the location of these quantum critical 
lines resulting from neglecting the effect of differential 
relaxation specific to the magnetic ordering and doping, 
and because the SGGA is not exact, is indicated by the 
broader yellow regions. 

The most notable feature of Fig. [I] is the asymmetry 
of electron and hole doping. The undoped material ex- 
hibits S-AFM ordering for all values of dpe-As (horizontal 
line, 5 = 0) and S-AFM ordering is stabilized by increas- 
ing dpe-As- On doping LaOFeAs with electrons (verti- 
cal line, d = do), the S-AFM ordering is weakened and 
a transition to NM ordering is observed at 5 ~ 0.18. 
The suppression of magnetic ordering by electron dop- 
ing is consistent with previous first-principles calcula- 
tions [25j [40] . We can use this phase diagram to explain 
the dependence of T c on external pressure. Experimen- 
tally it has been found that T c increases with external 
pressure for x < 0.14 and decreases for x > 0.14 for 
Smpi-^FjFeAs [15] . This behavior indicates the exis- 
tence of a QCP around x ~ 0.14. According to Fig.[l]the 
over-doped arsenide does not order magnetically while 
the under-doped (x < 0.14) material favours S-AFM or- 
der. External pressure will reduce dpe-As (white arrows) 
but its effect on T c depends on the doping. For under- 
doped SmfOi-^FjFeAs it destabilizes the S-AFM order- 
ing and pushes the system towards the QCL making it 
more susceptible to spin fluctuations. In agreement with 
observations [15 , we expect this to enhance T c . How- 
ever, the over-doped system is NM and pressure drives 
it away from the QCL leading to lower values of T c , as 
observed. When the structure is doped with holes, a tran- 
sition to C-AFM ordering occurs. We will discuss below 
how the asymmetry of hole and electron doping can be 
understood in terms of the density of states (DOS) of 
the parent compound close to the Fermi level. Our phase 
diagram makes a clear prediction for the pressure depen- 
dence of T c of hole doped iron arsenides, namely that it 
should be precisely the same as the electron-doped case 
electron-doped case. 

A second notable feature of Fig. [I] is the existence of a 
narrow region between NM and S-AFM states where the 



3 




Energy (eV) 



FIG. 2: Projected densities of states calculated for C-AFM 
ordered undoped and 25%-doped LaOFeAs; for S-AFM Fe or- 
dered undoped and 25%-doped LaOFeAs. The orbital char- 
acters are color-coded. The Fermi energy is shown by the 
vertical dotted line. 
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FIG. 3: Phase diagram of La[Oo.75Fo.25]FeAs as a function of 
lattice parameter a and inter-layer distance between Fe and 
As plane along c axis (^Fe-As)- The structure parameters 
obtained by geometry optimization for 25% doped Re com- 
pounds are plotted as yellow spheres. The color coding is the 
same as in Fig. [I] The white circle indicates the position of 
the tricritical point for a doping of 5 = 0.22. 



ground state is C-AFM for large dpe-As and high elec- 
tron doping. The line separating NM and C-AFM states 
meets the line separating C-AFM and S-AFM states for 
an electron doping of S ~ 22% and d^ e -As — 1.03do m a 
tricritical point (TCP) where NM, S-AFM, and C-AFM 
states coexist. For electron doping S > 22%, increasing 
^Fe-As generates a sequence of NM C-AFM -> S-AFM 
groundstates. To understand why this C-AFM region ex- 
ists, we examine the Fe projected DOS shown in Fig. [2] 
for C-AFM and S-AFM ordering for undoped and 25% 
electron doping. For S-AFM ordering, the Fermi energy 
in the undoped case is situated in a pseudogap for both 
spin channels. On doping with electrons, the pseudo- 
gap has to be crossed before a high density of (mainly 
minority-spin) states peak with d xy and d x i_ y i charac- 
ter can be populated. This is energetically unfavourable. 
For C-AFM ordering, no such pseudogap exists and the 
doping electrons can be accomodated in states close to 
the undoped Fermi level while the large value of c^Fe-As 
stabilizes the AFM states. It is the competition between 
these two mechanisms that results in the existence of a 
narrow region of C-AFM state. 

In the context of heavy fermion and cuprate supercon- 
ductivity, it has been suggested that spin fluctuations 
close to QCPs can give rise to the attractive interaction 
between carriers needed to mediate superconductivity 
[301 El [32]. The coexistence of NM, C-AFM, and S-AFM 
states at the TCP (S 22%;<ipe-As — 1.03c?o) suggests 
that spin fluctuations there may be stronger and more 
interesting. We therefore propose searching for higher 
values of T c in the vicinity of the TCP. This requires 
simultaneous achievement of large electron doping and 
negative pressures by suitable chemical substitutions. 

Chemical pressure can be exerted by replacing La with 



other rare earth atoms, the in-plane lattice constant a 
and unit cell volume decreasing monotonically with de- 
creasing Re size in the sequence La, Ce, Pr, Nd, Sm, Gd 
[33] . The effect of these substitutions has been studied 
experimentally resulting in the observation of a maxi- 
mum value of T c = 55K for SmOFeAs [7\ E E EED] • These 
substitutions change not only the lattice constants but 
also the internal parameter dp e -A s . To study how sen- 
sitively the QCPs depend on a, we construct a phase 
diagram in the (5 = 0.25, a, dpe-As) plane, Fig. [3] To 
achieve this high electron doping, we replace one out of 
four O atoms by a F atom in a y/2 x \[2 x 1 supercell 
in order to be able to more accurately describe distor- 
tions induced by the high concentration of the smaller F 
ion. The results agree well with those obtained on doping 
with the jellium model used to construct Fig. [I] 

The main features of this phase diagram are (i) that 
the effect of increasing a is similar to that of increas- 
ing <ipe-As : it stabilizes AFM ordering, and (ii) that the 
values of these structure parameters are not optimal. A 
small C-AFM "island" appears at the boundary between 
the NM and S-AFM "seas" in Fig. [3] consistent with 
Fig. [I] This C-AFM island represents an area where 
spin fluctuations are expected to be large. For a doping 
5 ~ 0.22 corresponding to the TCP, this island shrinks to 
a point located at d^ e -As = 1.23A, a = 4.04A, indicated 
in the figure by the white circle. We optimize the geome- 
tries of Re[O .75F .25]FeAs for Re =Ce, Pr, Nd, Pm, Sm, 
Gd, Tb, and Dy by energy minimization and plot the 
resulting structure parameters as yellow spheres in the 
figure [41]. Assuming that T c is proportional to proxim- 
ity to the TCP, we see that the Nd, Pm and Sm com- 
pounds should have the highest T c s in agreement with 
observation [7]. The experimental data for ReOi-^FeAs 
(Re=La, Ce, Pr, Nd, Sm) [8] have been interpreted as 
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evidence for a lattice contraction induced enhancement 
of T c . We note that replacing La with heavier Re ele- 
ments is accompanied not only by a reduction of a but 
also by an increase in dpe-As- From the phase diagram 
of Fig. [I] we attribute the higher T c primarily to the in- 
ecrease in dpe-As- We expect that the highest T c will 
occur for doped NdOFeAs and SmOFeAs and that re- 
placing Nd or Sm with heavier Re elements [10] will not 
lead to further improvement. 

The phase diagram of Fig. [l] suggests that the optimal 
electron doping for FeAs based superconductors is around 
22%. Electron doping can be reasonably well controlled 
experimentally by replacing O with F or by introducing 
oxygen vacancies. However Figs, [I] and [3] together indi- 
cate that currently achieved values of a and dpe-As are 
smaller than the optimal values and should be increased 
to achieve higher values of T c . One way of doing this 
would be to replace some As with the larger isovalent 
elements Sb or Bi or with Te. 

In conclusion, we propose that the electronic structures 
and magnetic properties of FeAs-based superconductors 
such as LaOFeAs, BaFe2As2 and LiFeAs are dominated 
by the effective doping and structural details of the FeAs 
layers. In a QCP scenario, the superconductivity is re- 
lated to spin fluctuations making it important to under- 
stand the magnetic phase diagram as a function of ex- 
perimentally accessible parameters. Currently available 
experimental data can be understood using the theoret- 
ically obtained phase diagrams we have presented giving 
us confidence that they can then be used to search for 
higher values of T c . 
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